Coronal rain as a marker for coronal heating mechanisms 
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ABSTRACT 

Reported observations in Ha, Ca II H and K or or other chromospheric lines 
of coronal rain trace back to the days of the Skylab mission. Offering a high 
contrast in intensity with respect to the background (either bright in emission if 
observed at the limb, or dark in absorption if observed on disk) these cool blobs 
are often observed falling down from high coronal heights above active regions. 
A physical explanation for this spectacular phenomenon has been put forward 
thanks to numerical simulations of loops with footpoint concentrated heating, a 
heating scenario in which cool condensations naturally form in the corona. This 
effect has been termed "catastrophic cooling" and is the predominant explana- 
tion for coronal rain. In this work we further investigate the link between this 
phenomenon and the heating mechanisms acting in the corona. We start by ana- 
lyzing observations of coronal rain at the limb in the Ca II H line performed by the 
Hinode satellite. We then compare the observations with 1.5-dimensional MHD 
simulations of loops being heated by small-scale discrete events concentrated to- 
wards the footpoints (that could come, for instance, from magnetic reconnection 
events), and by Alfven waves generated at the photosphere. It is found that if 
a loop is heated predominantly from Alfven waves coronal rain is inhibited due 
to the characteristic uniform heating they produce. Hence coronal rain may not 
only point to the spatial distribution of the heating in coronal loops but also to 
the agent of the heating itself. We thus propose coronal rain as a marker for 
coronal heating mechanisms. 
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1. Introduction 

Coronal loops are dynamical entities which exhibit heating and cooli ng; processes con- 



stant ly. Most of them are actually far from being in hydrostatic equilibrium (lAschwanden et al 



20011 ) . The dynamical nature of loops can manifest in observations in a variety of forms, 
among which propagating intensity variations is a common one. Indeed, intensity variations 
traveling along coronal loops have been frequently observed in many different wavelengths, 
in hot coronal EUV lines as well as chromospheric cool lines. The agents producing these 
intensity variations can be either propagating waves or flows along coronal loops, and ob- 
servers can have a hard time differentiating them despite their very different physical nature. 
Coronal rain is an example of intensity variations caused by flows. This spectacular phe- 
nomenon corresponds to cool plasma condensations (with chromospheric temperatures) in 
a hot environment (coronal temperatures) falling down along coronal loops, hence the very 
appropriate name. Coronal rain seems to be a common phenomenon above active regions, 
where loops are dense. Due to the low temperatures, the condensations appear as bright 
emissi on profiles in chromo spheric lines such as Ha or Ca II H and K when observed at the 



limb (IDe Groof et al.l 12004 . Fig. [T]) or as dark absorption profiles when observed on disk. 



If ob served in EUV wavelengths these propagating blobs appear as dark features (jSchrijver 



200l|). 



De Groof et al.l (12004 ); Ide Groof et al.l (120051 ) report simultaneous observations in EIT 



304 A, Ha images from Big Bear Solar Observatory, and in 171 A (TRACE) of intensity 
variations propagating along a coronal loop. Following a detailed analysis they put forward 
a series of points through which a differentiation between waves and flows (cool plasma con- 
densations, coronal rain) can be done. Propagating waves frequently reported in EIT/ SOHO 
195 A and TRACE 171 A appear as low min-to-max intensity variations as compared to 
coronal rain. Most observed waves correspond to sound modes propagating at constant speed 
corresponding to the local sound speed of the corona, ~ 150 km s . On the other hand, cool 
plasma condensations are seen to accelerate while falling along coronal loops to velocities 
above 100 km s _1 . Furthermore, waves have only been seen to propagate upwards, within 
the first 20 Mm of loops, after which they are usually damped, while cool condensations have 
only been seen to fall down along the loops from coronal heights. The falling speeds of the 
condensations have been reported to be lower than free fall speeds, due probably to the gas 
pressure along the loop, which increases considerably below the transition region. In some 
cases continuous flows from one footpoint to the other are also observed. 



Antiochos et al.l (119991 ) propose a common physical mechanism for coronal rain and 
prominence formation. Ranges for temperatures and densities seem to be shared by both 
mechanisms, as well as the location of formation, high in the corona. While a prominence 
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forms and is supported by the magnetic field for a time scale of days, coronal rain occurs on a 
time scale of minutes. It was shown that prominences may not need to rely on the geometry 
of the magnetic field to form (presence of "dips"), contrary to the general belief. They 
showed that a coronal loop whose heating is concentrated towards the footpoints is subject 
to a thermal instability in the corona, which dramatically cools down to chromospheric 
temperatures in time scales of minutes once the density and temperatures have reached 
critical values. This phenomenon was termed "catastrophic cooling" and has so far gained 
acceptance as possible explanation for coronal rain. This mechani sm occurs in a cyclic 
manner, matching reported periodicity of coronal rain observations (jSchrijverl l200ll ). This 
periodicity depends on geometrical aspects such as the loop length a nd the heating scal e 
height, and can be obtained even with a time independent heating (jMuller et al.l 120031 ) . 
As the heating is concentrated towards the footpoints a larger mass flux is input into the 
corona as compared to uniform heating. The density of the corona increases in time while 
the energy flux is kept constant due to the constant heating input at the footpoints. This 
causes the loop to reach a maximum temperature after which it starts decreasing due to 
the decreasing of the heating per unit mass. The lower temperatures and higher densities 
increase the radiative losses which further cool the corona. Eventually, in a time scale of 
hours, the loop reaches a critical state and a thermal instability sets in. Then, in a time 
scale of minutes the temperature and density respectively decrease and increase dramatically 
to chromospheric values. The cool condensation subsequently falls down subject to gravity 
and plasma pressure. The loop subsequently is depleted and reheats rapidly due to the low 
density. The cycle then repeats. We will refer to these cycles as "limit cycles" , as termed by 
flMuller et al"1l2003h . 



Observations seem to indicate that Coronal rain is a fairly common phenomenon of 
active regions, where loops are hot and dense. The coronal heating mechanism being a 
great unknown in solar physics, it is interesting to think about its underlying relationship 
to coronal rain. At a first glance coronal rain may appear as a random failure of the coro- 
nal heating mechanism to heat the loops in active regions. However, as previously stated, 
this phenomenon seems to act in the corona in a cyclic manner. Furthermore simulations 
have pointed out to the necessity of specific spa tially dependent heating in order to allow 
the catastrophic cooling leading to coronal rain (lAntiochos et al.lll999t iMuller et al.ll2003l ). 
This cooling phenomenon seems then to be deeply linked to the unknown coronal heating 
mechanism. Since it is a fairly easily observable phenomenon due to the large velocities and 
density variation (hence clear Doppler shifted emission or absorption profiles), it may act as 
a marker of the operating heating mechanism in the loop. It has been shown for instance that 
footpoint concentrated heating may lead to catast rophic cooling if th e heating scale height 
is sufficiently concentrated towards the footpoints (jMuller et al.ll2003l ). Uniform heating on 
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the other hand fails to reproduce the phenomenon since the heating rate per unit mass needs 
to decrease in time locally in the corona in order to allow the thermal instability to set in. 
It was further shown that catastrophic cooling does not need time dependent heating. In 
other words, it can happen even with a constant heating function. The footpoint concen- 
trated heating function to which simulations point to matches the observation al evidence of 



coron al loops above active regions for being mainly heated at their footpoints (lAschwanden 



20011 ) . which sets most loops out of hydrostatic equilibrium. Further evidence of this fact 
has been found by lHara et al.l (120081 ) using the Hinode/EIS instrument, which shows that 
active region loops exhibit upflow motions and enhanced nonthermal velocities. Possible 
unresolved high-speed upflows were also found, fitting in the footpoint concentrated heating 
scenario. 

Many heating mechanisms have been proposed as candidates for heating the solar corona 
up to the observed few million degree temperatures. In this context a large emphasis has 
been put on the search for Alfven waves in the solar corona. Theoretically, they can be easily 
generated in the photosphere by the co nstant turbulent convective motions, whi ch inputs 
large amounts of energy into the waves (iMuller et al.lll994t IChoudhuri et al.lll993l ). Having 
magnetic tension as its restoring force the Alfven waves can travel less affected by the large 
transition region gradients with respect to other modes. Also, when traveling along thin 
magnetic flux tubes they are cut-off free since they are not coupled to gravity (Musielak et 
al. 2007)0. Alfven waves generated in the photosphere are thus able to carry sufficient en- 
ergy into the corona to compensate the losses due to radiation and conduction, and, if given 



a suitable dissipation mechanism, heat the p 



peratures (lUchida fc Kaburakilll974l; IWentzel 



Ruderman et al 



199 



7; 



Kudoh &; Shibata 



asma to the high mi 



1974 



Holl weg et al. 



lion degree coronal tem- 



1982; IPoedts et alJll989 



1999; lAntolin fc Shibatal 120091 ) and power the solar 



wind ( jSuzuki fc Inutsukal 120061 ; ICranmer et al.l 120071 ) . The main problem faced by Alfven 
wave heating is actually to find a suitable dissipation mechanism. Being an incompress- 
ible wave it must rely on a mechanism by which to convert the magnetic energy into heat. 
Several dissipation mechan isms have been pro posed, such as parametric decay ([Goldstein 



1978; Terasawa et a 



19861 ). mode c onversion dHollweg et al 



Moriyasu et al.ll2004l ). p hase mixing dHewaerts fc Priest 



or re sonant absorption (llonsonlll978l ; iHollweg 



1984 



1983 



1982 



Kudoh fc Shibata 



lOfman fc Aschwanden 



1999 



2002 



Poedts et al.l 11989c lErdelyi fc Goossens 



19951 ). The main difficulty lies in dissipating sufficient amounts of energy in the correct 
time and space scales. For more discussion regarding this issue the reader can consult for 



Werth et al.l (|2009. in preparation! ) have pointed out however that the assertion made bv lMusielak et al. 
(J2007]) is valid only when the temperature in the flux tube does not differ from that of the external plasma. 
When this is not the case a cut-off frequency is introduced. 
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instance iKlimchukl (120061 ). lErdelyi fc Ballail (120071 ) and lAschwandenl (120041 ) . Works con- 
sidering Alfven wave heating as coronal heating mechanism have shown that the obtained 
coronae are uniformly heated (IMoriyasu et al.l 120041 ; lAntolin et al.l 120081 ; lAntolin fc Shibata 
20091 ; ISuzuki Inutsukall2006l ). In these works the heating issues from shocks of longitudinal 
modes (mainly slow modes) from mode conversion of the Alfven waves due to the density 
fluctuation, wave-to-wave interaction and deformation of the wave shape during propaga- 
tion. The coronal loops issuing from Alfven wave heating are found to satisfy quite well the 
RTV scaling law ( Rosner et al.l [l978l ) which quantifies the heating uniformity in the loops. 
This result would point then towards an inhibition of coronal rain if Alfven wave heating is 
predominant in the loop. It is this idea that is addressed in this work. 

Another promising coronal heating candidate mechanism is the na noflare reconn ection 



heating model. The nanoflare reconnection process was first suggested by lParkerl (119881 ). who 
considered a magnetic flux tube as being composed by a myriad of magnetic field lines braided 
into each other by continuous footpoint shuffling. Many current sheets in the magnetic 
flux tube would be created randomly along the tube that would lead to many magnetic 
reconnection events, releasing energy impulsively and sporadically in small quantities of the 
order of 10 24 erg or less (nanoflares) . Parker's original idea was a nanoflare reconnection 
heating acting uniformly in the corona but it was later proposed to be concentrated towards 



( Klimchuk 


2006; 


Aschwanden 


2001) 



to be generated, and the energy imparted into Alfven waves is a matter of debate. The 
impart ed energy ma y well depend on the location in the atmosphere of the reconnection 



event. IParkerl (119911 ) suggested a model in which 20 % of the energ y released by reco nnection 
events in the solar corona is transfered as a form of Alfven wave ] Yokoyamal (119981 ) studied 
the problem simulating reconnection in the corona, and found that less than 10 % of the total 
released energy goes into Alfven waves. This result i s simil ar to the 2-D simulation results 
of photospheric reconnection by iTakeuchi fc Shibatal (120011 ). in which it is shown that the 
energy flux carried by the slow magnetoacoustic waves is one order of magnitude higher that 
the energy flux carried by Alfven waves. On the other hand, recent simulations by Kigure et 
al. (private communication) show that the fraction of Alfven wave energy flux in the total 
released magnetic energy during reconnection in low (3 plasmas may be significant (more 
than 50 %). Since the observed ubiquitous inten sity bursts (n anoflares) are thought to play 
an important role in the heating of the corona ([Hudson! Il99ll ) and since they are generally 
assumed to be a signature of magnetic reconnection it is then crucial to deter mine the energy 



going into the Alfven waves during the reconnection process. Moreover, IMoriyasu et al. 



(120041 ) has shown that the observed spiky intensity profiles due to impulsive releases of energy 
may actually be a signature of Alfven waves. It was found that due to nonlinear effects Alfven 
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waves can convert into slow and fast magnetoacoustic modes which then steepen into shocks 
and heat the plasma to coronal temperatures balancing losses due to thermal conduction 
and radiation. The shock heating due to the conversion of Alfven waves was found to be 
episodic, impulsive and uniformly distributed thr oughout the corona, p roducing an X-ray 
intensity profile that matches observations. Hence, iMoriyasu et al.l (120041 ) proposed that the 
observed nanoflares may not be directly related to reconnection but rather to Alfven waves. 

Differentiating Alfven wave heating from nanoflare reconnection heating during observa- 
tions is one of the main task s need e d in order to so l ve the coronal heating problem. Following 
the work of IMoriyasu et al.l (120041 ). lAntolin et al.l (120081 ) have compared both heating mech- 
anisms by studying the hydrodynamic response of a loop subject to both kinds of heating. It 
was found that Alfven waves lead to a dynamic, uniformly heated corona with steep power 
law indexes (issuing from statistics of heating events) while nanoflare-reconnection heating 
leads to lower dynamics (besides the times when catastrophic cooling takes place in the case 
of footpoint concentrated heating) and shallow power laws. It was further found that foot- 
point nanoflare heating (namely, nanoflare reconnection heating concentrated towards the 
footpoints of loops) leads to hot upflows (as observed in the FeXV 284.16 A line) due to 
the plasma being heated rapidly towards the footpoints before being ejected into the corona, 
while Alfven wave heating leads to hot downflows due to the plasma achieving the maximum 
temperatures in the corona (rather than at the footpoints) a n d bei ng carried back to the 
footpoints by the strong shocks generated there lAntolin et al.l (120091 ) . In this work we pro- 
pose coronal rain as another observational signature through which both heating mechanisms 
can be distinguished. 

We start first by reporting limb observations of coronal rain from Hinode /SOT in the 
Call H line. The velocities and shapes of the falling condensations are analyzed. With a 
1.5-dimensional code we then proceed to model a coronal loop being subject to a heating 
mechanism that is concentrated towards the footpoints, such as the nanoflare reconnection 



heating model (as proposed by lAschwandenl 120011 ; iKlimchukl 120061 ) . In the case considered 
catastrophic cooling happens 3 times and we select the first event to analyze and compare 
with our observations. Next we proceed by generating Alfven waves at the photosphere and 
gradually increase the amplitude. As the heating from the waves becomes non-negligible with 
respect to the heating flux from the nanoflare heating events catastrophic cooling is inhibited 
and the loop reaches a thermal equilibrium state. This result has important consequences 
since it points to the conclusion that Alfven waves are a non predominant heating mechanism 
in active region loops. 

The work is organized as follows. In £j2]we report observations of coronal rain in the Ca 
II H line performed with Hinode/SOT. In §3] we introduce the 1.5-dimensional MHD model 
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in which our loop is based on and discuss the heating models of the loop. In §Hwe present 
the results of footpoint heating and analyze a typical case of catastrophic cooling. The effect 
of Alfven waves on the thermal stability of the loop is also studied. In §5] we discuss the 
results in the context of coronal heating and conclude the work. 



2. Observations of coronal rain with Hinode/SOT 

In this section we report high resolution observations of coronal rain at the limb per- 
formed by the Solar Optical Telescope on board of Hinode in a 0.3-nm broadband region 
centered at 396.8 nm, the H-line spectral feature of singly ionized calcium (Ca II H line). 
The observation was performed on 9 November 2006 from 19:33 to 20:44 UT with a cadence 
of 15 s, and focused on NOAA AR 10921 on the west solar limb. The corresponding field of 
view is 80 Mm x 40 Mm. The Ca II H line is a chromospheric line typically showing plasmas 
with temperatures on the order of 20,000 K. In the data set, apart from the active region on 
disk, many interesting structures can be seen over the lim b, such as spicu l es, pr ominences 



and coronal rain (Fig. [T]). The same data set was used by lOkamoto et al.l (120071 ) to report 
observations of transverse magnetohydrodynamic waves propagating in the observed promi- 
nences. The later show complex horizontal threads displaying continuous horizontal motions, 
and appear to be located on the background of the images. Coronal loops exhibiting coronal 
rain appear to be located on the foreground of the prominences. The cool condensations 
falling down the loops form close to the apex at heights of ~ 40 Mm above the surface, and 
do not seem to be linked to the prominences. 

The condensations conforming coronal rain trace the geometry of the coronal loops and 
appear to accelerate as they fall down to the chromosphere under the action of gravity. By 
fitting the loops with a curve we can trace the condensations along their way down and 
calculate the corresponding velocities. In Fig. [2] we show length-time diagrams of 3 loops 
exhibiting coronal rain located on the upper left side of Fig. [TJ The estimated velocities 
on the top part of the loops (near the apex), and close to the footpoints are displayed in 
Table ??. We also calculate the free fall velocity that the condensations would have under 
the action of gravity alone and no initial velocity (at the top part of the loops). We can see 
that close to the top parts of the loop the condensations have velocities around 50 - 70 km 
s _1 , from which we estimate the apexes of the loops to be at 5 - 10 Mm above the top of the 
figures (assuming that the condensations start with zero velocity at the apex). The lengths of 
these loops should then be close to 100 Mm. At the footpoints of the loops the condensations 
have their highest velocities, roughly between 115 and 140 km s _1 . The calculated free fall 
velocities seem to match these values, however, they have been calculated assuming no initial 
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velocity. We can then expect a discrepancy on the order of the velocity at the top (~ 50 km 
s -1 ). Hence, the condensations fall at velocities which are lower than the free fall velo cities, 
as has been reported previously for coronal rain observations ( IDe Groof et al.ll2004l ). The 
decelerations may be caused by internal pressure changes in the loops, as suggested by the 
results of the simulations reported in section HI For the loop in the right panel we have a 
clear deceleration of the coronal rain from ~ 70 to ~ 15 km s _1 halfway down along the loop. 
It then seems to accelerate again to a speed above 100 km s _1 , higher than the corresponding 
free fall velocity. For the loop in the middle panel we can see an upward motion of coronal 
rain at ~ 25 km s _1 close to the top part of the loop. This may be caused by a cooling 
and heating process acting upward in that part of the loop, thus resulting in an apparent 
upward motion of the coronal rain. However, it may also be a true motion caused by a 
pressure change inside the loop at that height, such as from propagating acoustic shocks, as 
the simulations performed in H] suggest. 



3. Simulation setup 



In order to simulate coronal rain we follow the model of lAntiochos et al.l (119991 ). in 
which it is shown that cool condensations can dynamically form in the corona resulting from 
footpoint heating of the loops. This mechanism is known as catastrophic cooling and is 
further explained in this section. 



3.1. Model 



We consider a magnetic flux tube (loop) of 100 Mm in length, roughly the same length 
as the loops exhibiting coronal rain in the observations reported here. The geometry of the 
loop takes into account the area, which considers the predicted expansion of magnetic flux 
in the photosphere and chromosphere, displaying an area ratio between the corona and the 
photosphere of 1000. As discussed in the introduction catastrophic cooling as proposed by 
Antiochos et al.l (Il999l ) needs the heating to be concentrated towards the footpoints of loops. 
In the case of footpoint heating without the generation of Alfven waves in the photosphere 
the plasma motion is governed by the usual ID HD equations for the conservation of mass, 
momentum and energy. The model in this case is the same as the heating mo del considered 



for na noflare reconnection with heating concentrated towards the footpoints in lAntolin et al. 



(120081 ). and is further discussed below. When Alfven waves are considered the model gains 
the azimuthal comp onent and is the same model as the model for Alfven wave heating in 
Antolin et all (l2008h . 
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We write the ID HD equations for the conservation of mass, momentum and energy in 
the following way: 



the mass conservation equation: 



the momentum equation: 



and the energy equation: 



dp dp n 3 /» 

m +v dS = - pB d:s Kb 



dv ^ dv 1 dp _ 
dt ds pds 



de de , ^ „ d ( v \ R-S-H 1 d ( 2 dT\ 



d-t +v dS = ~ {l - 1)eB dS (b) p + pV^dS [ r K ^J ; (3) 

where 

p = p—T, e= --. 4 

m 7 — 1 p 

In the above equations CQ)-© s measures the distance along the flux tube (central field line) 
and r is the radius of the tube, p, p, v and e are, respectively, density, pressure, velocity 
along the loop and internal energy; B is the magnetic field along the loop and is a function 
of r alone, B = B (ro/r) 2 , where B is the value of the magnetic field at the photosphere 
and r = 200 km is the initial radius of the loop; hs is the Boltzmann constant and 7 is the 
ratio of specific heats for a monatomic gas, taken to be 5/3. The g s is the effective gravity 
along the loop and is given by 

9s = 9& cos (|-7rJ , (5) 

where g & = 2.74 x 10 4 cm s -2 is the gravity at the base and L is the total length of the loop. 

We assume an inviscid perfectly conducting fully ionized plasma. The effects of thermal 
conduction and radiation are taken into account, where the Spitzer conductivity correspond- 
ing to a fully ionized plasma is considered, and radiative losses are defined as 

R(T)=n e n p Q(T) = —Q(T). (6) 

Here, n = n e + n p is the total particle number density (n e and n p are, respectively, the 
electron and proton number densities, and we assume n e = n p = p/m to satisfy plasma 
neutralit y, with m the proton mass) and Q( T) is the radiative loss function for optically thin 



plasmas (jLandini &: Monsignori Fossilll990l ) which is approximated with analyt i cal fu nctions 



of the form Q(T) = xT 1 ■ We take the same approximation as in iHori et al.l (119971 . please 



refer to their Table 1). Apart from the optically thick regions in the low solar atmosphere, it 
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is possible to have optically thick plasma in some frequencies high up in the corona by means 
of catastrophic cooling of loops which would lead to coronal rain. The dynamical state of the 
plasma in this situation (catastrophic cooling events happen in the time scale of minutes) 
would produce non-equilibrium ionization effects which should be taken into account by 
solving the ionization rate equations in a self-consistent way (by coupling the equations to 
the hydrodynamic equations through the energy equation). Here we are then assuming that 
the radiation fields in all directions and all frequencies and the level populations do not 
affect the ionization level of the plasma. Since coronal rain disappears in the same time scale 
(minutes), we may expect the plasma that becomes optically thicker not to affect considerably 
the energy equation. We consider however an approximation to optically thick radiation. For 
temperatures below 4 x 10 4 K we assume that the plasma inside the loop becomes optically 
thick. In this case , the radiative losses R can be approximated by R(p) = 4.9 x 10 9 p 



( lAnderson fc Athayi Il989l ). In equation §3§ the heating term S has a constant non-zero 
value which is non-negligible only when the atmosphere becomes optically thick. Its purpose 
is mainly for maintaining the initial temperature distribution of the loop. Here TC denotes 
the heating function in the loop, which corresponds to a nanoflare heating model and is 
presented in the next section. 



For generating Alfven waves in the loop we follow the same model as in lAntolin et al. 



( 120081 ). Random torque motions are produced in the photosphere, which generates Alfven 
waves with a white noise spectrum in frequency. We adopt this model instead of a monochro- 
matic wave generator since we consider that the buffeting of magnetic field lines by convective 
motions has a turbulent nature, thus le a ding t o random motions. For further details about 



this model please refer to lAntolin et al.l (120081 ). 



3.2. Nanoflare heating function 



As shown by lAntiochos et al.l (119991 ). in order for catastrophic cooling to happen we have 
to apply a heating mechanism that is concentrated towards the footpoints of the loop. There 
may be many proposed heating mechanisms which can act preferentially towards the foot- 
points of loops. Here we will assume that the loop is subject to "fo o tpoin t nanoflare" heating 
from the nanoflare reconnection model described in lAntolin et al.l (120081 ) . In this picture we 
assume that the energy imparted into Alfven waves from reconnection events is low and can 
be neglected relative to the imparted energy to the slow modes. Hence, in this picture the 
corona would be heated mainly by the accumulation of numerous nanoflares coming from re- 
connection events and by slow magnetoacousti c shocks. Hydrodynamic modeling of nanoflare 



heating has already been done in the past (IWalsh et al.l Il997l ; ICargill fc Klimchuki 12004 



Patsourakos fe Klimchukll2005l ; iTaroyan et al. 2006; Mcndoza- Bricefio et al. I l2005h . The nanoflare 



model considered here is similar to the model of ITaroyan et al.l (120061 ) with respect to the 
heating function TC in equation In the present case we assume that heating events sim- 
ulating reconnection events (leading to nanoflares) occur towards the footpoints of the loop. 
These are input randomly as artificial perturbations in the internal energy of the gas (thus 
generating only slow modes). The heating rate due to the nanoflares is represented as 

n 

n = ^2?u{t,8) (7) 

i=l 

where TCi(t,s), i=l,...,n are the discrete episodic heating events, and n is the total number 
of events. 

£ sin(^)exp(-^), tl <t<U + n; 
0, otherwise, 



Hi(t, s) 



where Eq is the maximum volumetric heating and Sh is the heating scale length. The offset 
time ti, the maximum duration Tj, and the location Sj of each event are randomly distributed 
in the following ranges: 

ti ^ [0) ^totaz]; Ti £ [0; T~max\i &i ^ [^minj ^mra] l s _J[^" / ^maxi L Smin]) (9) 

where t tota i is the total simulation time and s min (s max ) define the lower (upper) boundaries 
of the range in the loop where heating events occur. 

In order to set the values to the parameters of the heating function, equations ([Tj)-©, an 
estimate of the nanoflare duration time is needed. One of the hardest parameters to estimate 
in magnetic reconnection theory is the thickness of the current sheet, i.e. the length across 
the reconnection region. If this parameter is of the order of ~ 1000 km, the timescale of a 
(small) reconnection event leading to a nanoflare should oscillate between 1 and 10 s, since 
the order of the Alfven speed in the chromosphere and in the corona is, respectively, ~ 100 
and ~ 1000 km s _1 . This value, however, is not established. Different values have been 
tried for the parameters of the heating function defined in equations (J7]) - ©• Since the 
purpose of this work is not to study the ranges in which catastrophic cooling happens we 
will limit ourselves in the present model to present a typical case in which it happens. For 
this case we have the maximum duration time of a heating event r max = 40 s, the heating 
scale length Sh = 1000 km, the maximum volumetric heating Eq = 0.5 erg cm -3 s -1 , a 
frequency of heating events of 1 each 50 s, the upper and lower boundaries of the ranges in 
which heating occurs {s m i n = 2, s max = 20} Mm. These parameters set a mean energy per 
event of 1.9 x 10 26 erg and a mean energy flux of 2.5 x 10 7 erg cm -2 s _1 . 
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3.3. Initial conditions and numerical code 



The loop is assumed to follow hydrostatic pressure balance in the subphotospheric region 
and in the photosphere up to a height of 4Hq = 800 km, where Hq is the pressure scale height 
at z — 0. For the rest of the loop, density decreases as p oc h~ 4 , where h is the height from 
the base of the loop. This is based on the work by IShibata et al. rll989al lbh. in which the 
results of 2D MHD simulations of emerging flux by Parker instability exhibit such pressure 
distribution. The initial temperature all along the loop is set at T = 10 4 K. The density 
at the photosphere (z = 0) is set at po = 2.53 x 10~ 7 g cm -3 , and, correspondingly, the 
photospheric pressure is po = 2.09 x 10 5 dyn cm~ 2 . The plasma (3 parameter is chosen to be 
unity at z = 0, setting the photospheric magnetic to B Q = 2.29 x 10 3 G. The value of the 
magnetic field at the top of the loop is then B top = 2.29 G. 

The spatial resolution in the numerical scheme is set to 5 km up to a height of ~16000 
km. Then, the grid size slowly increases until it reaches a size of 20 km in the corona. 
The size is then kept constant up to the apex of the loop. The total grid number is 10000. 
We take rigid wall boundary conditions at the photos phere. The numerica l scheme adopted 
is the CIP ( c ubic interpolated propagation) scheme (jYabe fc Aokil Il99lh . Please refer to 
Kudoh et al.l (119981 ) for details about the application of these scheme. The total time of the 
simulation is 568 minutes. 



4. Simulation results 

We first perform the simulation of the loop being heated only by the events from the 
"footpoint nanoflare" model, that is, without Alfven waves. We then allow the generation of 
Alfven waves at the footpoints of the loop and analyze the effect on the catastrophic cooling 
events. 



4.1. Footpoint nanoflare heating 

Due to the large energy flux from the heating events the corona is formed rapidly in 
the considered footpoint nanoflare model (in about 20 min). The mean temperature in the 
corona over the entire simulation time is (T) ~ 1.4 x 10 6 K, with a maximum temperature 
of (T) ~ 3.9 x 10 6 K. The mean density in the corona is high, (n) ~ 1.3 x 10 9 cm -3 , 
characteristic of a dense active region loop. As the heating events occur randomly in time 
we have a uniform heating input into the loop. Since the heating events occur close to the 
footpoints, chromospheric matter is constantly being push upwards into the loop, increasing 
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the density in the corona. Fig. [3] is a phase diagram of the mean temperature and the 
corresponding mean density in the corona in time. Arrows indicate the time direction and 
curve styles indicate different cycles the loop experiences. In the present case 3 cycles can 
be noticed, each one lasting roughly 170 min. The dotted curve in Fig. [3] corresponds to 
the initial phase of the simulation, then the 3 cycles start, denoted by solid, dashed and 
dot-dashed curves (blue, green and red curves in the online version), respectively in time. 
A cycle is composed of 4 distinctive phases. First, a phase in which the temperature of 
the corona increases rapidly and the density is roughly constant. Then follows a phase of 
constant temperature and slow density increase. In the third phase the temperature in the 
corona slowly decreases and the density is roughly constant. The last phase is marked by a 
dramatic decrease of temperature which can happen either locally in the corona or globally 
(entire collapse of corona), accompanied by a dramatic increase of den sity (at one o r more 



locations in the corona). These cycles have been termed "limit cycles" by lMuller et al.l (120031 ) 
and can be understood as follows. 

Initially, the density in the corona is low, since gravity depletes the loop before having a 
considerable injection of mass from the heating mechanism. The corona is thus easily heated 
to high temperatures from the footpoint heating events (phase 1). These events continuously 
inject material into the corona, thus slowly increasing its density. Since the energy flux to 
the corona is kept constant, the high temperatures can only be kept for a specific density 
range (phase 2). Then the mean temperature starts decreasing due to the steady decreasing 
heating rate per unit mass (phase 3). The lower temperature also increases the (optically 
thin) radiative losses, accelerating the cooling of the corona. The loop then reaches a critical 
state in which its density is too high and the temperature too low. A thermal instability 
follows due to the large radiation increase at low temperatures (becoming optically thick) 
in just the same way the transition region forms (phase 4). The density (temperature) then 
rapidly increases (decreases) to chromospheric values in a time scale of minutes, thus leading 
to a catastrophic cooling event. The thermal instability can happen locally in the corona, 
thus forming a dense and cool blob which subsequently falls down due to gravity, or can have 
a more global character, case in which the entire corona collapses and several dense blobs 
are formed. The loop is then evacuated and gets rapidly reheated due to the low density 
and the constant heating input. The cycle thus restarts. In the cycle denoted by the solid 
line (blue line in the online version) of Fig. [3] the catastrophic cooling occurs locally in the 
corona, while in the subsequent cycles it occurs globally. As soon as the condensation forms, 
its density increases rapidly and it grows larger. This is due to the constant heating events 
happening towards the footpoints constantly pushing plasma upwards which accumulates in 
the condensation. This can be seen in Fig. HJ where the evolution in time of the density along 
the loop is shown. The (acoustic) shocks created by the heating events can be followed from 
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the transition region. Some of them also form small condensations while propagating, before 
colliding with the bigger original condensation. It can be noticed that the blob experiences 
a change of direction at the loop top, going first towards the left footpoint and then towards 
the right footpoint. A shock collides with the blob and is reflected just at the time in which 
the blob changes direction. Another factor of the change of direction of the blob is the lower 
gas pressure region on the right side of the blob, as compared to the left side. Hence, these 
condensations are not only subject to gravity but also to the local changes of gas pressure 
in the loops. This mechanism may explain the observed deceleration motion and upwards 
motion of coronal rain in Fig. [2j When the blob falls down to the chromosphere it experiences 
a strong deceleration by the higher density region, the transition region is left oscillating up 
and down a couple of times. 

In the upper panel of Fig. we track down the condensation as it falls down to the 
chromosphere and plot its velocity along the loop in time. The initial motion towards the 
left footpoint has a maximum speed of almost ~ 30 km s _1 before changing direction and 
accelerating towards the right footpoint under the action of gravity. The maximum speed is 
almost 120 km s _1 close to the footpoint, before being decelerated by the high gas pressure of 
the chromosphere. The obtained values for the velocity match well the velocities of coronal 
rain we have found from the observations with Hinode/ SOT. The lower panel shows the 
width of the condensation as it falls down the loop. We can see a general tendency to 
elongate, as some parts of the blob may accelerate faster due to the variation of the effective 
gravity along the loop, and also to the large accumulation of matter in the condensation. 
This elongation is observed for the loops exhi biting coronal r ain in the observations reported 



here, and has also been previously reported (jSchrijverl 120011 ). 



4.2. Footpoint nanoflare heating and Alfven waves 



In order to analyze correctly the effect of Alfven waves on the catastrophic cooling events 
we first consider the case in which we only have Alfven waves and no nanoflare reconnection 
heating. Fig.[6]shows the phase diagram of mean temperature and mean density of the corona 
in this case. We can see that as the simulation evolves the mean temperature and density 
in the corona converge rapidly to roughly constant values, seen as an attractor in the phase 
diagram. Indeed, after one fifth of the total simulation tim e (142 min) the temp eratures and 
densities of the corona stay roughly constant. As seen in lAntolin et al.l (120081 ) the corona 
that issues from Alfven wave heating is uniform and steady, and satisfies the RTV scaling 
law. 



We then consider an hybrid model in which both heating mechanisms are present. 
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Heating events simulating reconnection events happen close to the footpoints, and Alfven 
waves with a white noise spectrum are generated in the photosphere. We further allow the 
amplitude of the waves to increase in time, as shown in Fig. [7J At the beginning of the 
simulation the energy flux from the waves is negligible with respect to the energy flux from 
the nanoflare reconnection events . Indeed, Alfy e n way es have an amplitude lower than 0.3 
km s _1 , which from the study in lAntolin et al.l (120081 ) (see Fig. 4 in that paper) we know 
it's not enough to produce a hot corona. In the second half of the simulation the waves 
have an amplitude larger than 1 km s _1 (reaching ~ 2 km s" 1 by the end of the simulation), 
which, according to Fig. 4 in that paper, is enough to produce a hot corona. In Fig. [8] we 
plot the corresponding phase diagram for this case. We can see that as the amplitude of the 
waves increases the limit cycles get smaller and disappear. Mean temperatures and densities 
finally converge to roughly constant values. Hence the loop becomes uniformly heated as 
the heating from the Alfven waves is no longer negligible. 



Discussion and conclusions 



Observations in chromospheric lines such as Ha or Ca II H and K seem to show that 



coron al rain is a phenomenon exclusively of active regions (jSchrijverl l200ll ; |Pe Groof et al. 



200 41), where loops are dense and heating appears to be conc entrated towards the footpoints 



( lAschwanden et al.ll200ll ; lAschwandenll2001 



Hara et al.ll2008l ). In this work we have reported 



observations with the Hinode/SOT telescope in the Ca II H line of coronal rain occurring 
over an active region, and compared with results of simulations of loops that undergo catas- 
trophic cooling. We have found that catastrophic cooling of coronal loops satisfactorily 
reproduces the main observed features of coronal rain. The loops are prefere ntially heated 
towa rds the footpoints and are subject to cycles ("limit cycles", as termed by iMiiller et al. 



20031 ) in which they rapidly cool down and then reheat, they get dense and then deplete. 
The constant heating input at the footpoints of the loops produces coronae out of hydro- 
static equilibrium. The coronal density increases in time causing a gradual decrease of the 
temperature and increase of the radiative losses. When the temperatures are sufficiently low 
radiative losses are dramatically increased since the condensation becomes optically thick 
and radiates considerably more. Catastrophic cooling sets in, either locally in the corona or 
globally, case in which the entire corona is cooled down to chromospheric values. Dense con- 
densations of cool plasma form at coronal heights, which subsequently fall down by gravity. 
The loop is then depleted and the cycle starts again. 



The characteristics of our condensation match well the main characteristics of the coro- 
nal rain from the reported observations. The velocities are in the same range as would be 
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expected from coronal rain forming in roughly the same locations of loops having essentially 
the same length. We obtain an elongation of the condensation in the simulation as it gets 
denser and the effective gravity increases along its way down the loop. This elongation is 
observed as wel l in the present Hinode observations and has previously also been reported 
Schrijverl (120011 ). The temperatures and densities of the resulting condensation have chro- 
mospheric values, which suggest that they would emit radiation in lines such as Ha or Ca 
II H or K, as in the observations with Hinode/SOT or the Swedish Solar Telescope. In a 
future paper we will test this idea by synthesizing the Ha line profile based on self-consistent 
solution of the ionization rate equations coupled with the hydrodynamic equations. 

As suggested by the simulations coronal rain is subject not only to gravity but also 
to the pressure changes inside the loop. Catastrophic cooling is the abrupt loss of thermal 
equilibrium in the corona, through which large pressure changes along the loop can occur 
and which can drive strong flows and shocks along the loop. In the catastrophic cooling 
event from the simulation reported here an upward motion of the condensation is obtained 
due to the local changes in pressure from the flows and the shocks. This may be a possible 
explanation to the upward motion of coronal rain as reported from the observations with 
Hinode/SOT in the Ca II H line. 

If coronal rain is indeed the consequence of the catastrophic cooling mechanism we then 
must have a heating mechanism acting preferentially towards the footpoints in loops where 
coronal rain is observed, i.e. active region loops. We have seen that Alfven waves generated at 
the footpoints of loops produce uniform coronae and are thus unable to reproduce phenomena 
such as coronal rain. Furthermore, when Alfven waves are present in a loop and have enough 
energy flux to heat and maintain a hot corona the catastrophic cooling events are inhibited. 
The loop then converges to a uniform and steady state. Coronal loops in active regions show 
a recurrent occurrence of coronal rain. They are dynamical entities showing heating and 
cooling processes at all times. Our results indicate that having Alfven wave heating as the 
main heating mechanism in loops means the absence of coronal rain. These results indicate 
then that Alfven wave heating may not be the principal heating mechanism for coronal loops 
in active regions. 



Hara et al.l (120081 ). using the Hinode/EIS instrument, have found upflow motions and 



enhanced nonthermal velocities in the hot lines of Fe XIV 274 and Fe XV 284 in active 



region loop s. Possible u n resolv ed high-speed upflows were also found. In lAntolin et al. 



( 120081 ) and lAntolin et al.l (120091 ) we found that footpoint or uniform heating coming from 
nanoflare reconnection exhibit hot upflows, thus fitting in the observational scenario of active 
regions, while Alfven wave heating was found to exhibit hot downflows, which may fit in the 
observational scenario of quiet Sun regions (jChae et al.lll998l ; iBrosius et al.l 120071 ). further 
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supporting the present conclusions. In lAntolin fc Shibatal (120091 ) we have found that Alfven 



wave heating is effective only in thick loops (with area expansions between photosphere and 
corona higher than 600) and long loops (with lengths of the corona above 50 Mm), a scenario 
which may not fit in active regions, where loops exhibit low area expansions due to the high 
magnetic field filling factors in those regions. Hence Alfven waves may play an important 
role in the heating of Quiet Sun regions (rather than active regions), where loops are often 
long, expand more than in active regions, kG (or higher) bright points are ubiquitous and 
coronal rain appears to be absent. 
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Fig. 1. — High resolution Hinode/SOT image at the limb in the Ca II H line. The entire 
field of view corresponds to roughly 80 Mm x 40 Mm. The observation was performed 
on November 9, 2006, between 19:33 to 20:44 UT with a cadence of 15 s. In order to be 
able to see the fainter corona the intensity of the photosphere (disk) has been decreased. 
Also, the intensity of the limb structures (spicules) is saturated. The sunspot visible on disk 
corresponds to the main sunspot of NOAA AR 10921. Above the limb, cloud-like prominence 
structures and coronal rain can be observed. The intensity of coronal rain in the Ca II H 
line is about 1% of the on-disk photospheric intensity. 
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Fig. 2. — Length-time diagrams for 3 loops exhibiting coronal rain on the upper left of Fig. 
CD The loops are traced by quadratic interpolation. The y axis denotes the length along 
the loops. Coronal rain, which appears brighter in the Ca II H line intensity, generally 
accelerates downward under the action of gravity. However, deceleration and even upward 
motions can also be observed (right and middle panels respectively). 
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Fig. 3. — Phase diagram of mean temperature and mean density of the corona. Arrows show 
the time direction and solid, dashed and dot-dashed curves denote the limit cycles (blue, 
green and red in the online version). The circle, triangle, square and lozenge denote the end 
of these cycles, respectively. The dotted curve corresponds to the start of the simulation. 
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Fig. 4. — Density map along the loop. Catastrophic cooling occurs forming a cool and dense 
condensation close to the apex of the loop, which falls down with increasing speed due to 
gravity. The motion of the condensation is also subject to the local gas pressure which varies 
considerably due to the acoustic shocks produced by the heating events at the footpoints 
of the loop. The traces of the propagating shocks are clearly observed. Note the strong 
deceleration of the blob as it enters the chromosphere, and the following oscillation of the 
transition region. 
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Fig. 5. — Upper panel: velocity of the condensation along the loop with respect to time. 
t=0 min corresponds to the formation time of the condensation. Lower panel: Width of the 
condensation with respect to time since its formation time in the corona. 
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Fig. 6. — Phase diagram of the mean temperature and density in the corona in the case 
of a loop heated by Alfven waves. The arrow indicates the time direction. The times 
corresponding to the circle and the triangle (end of simulation) are indicated. Limit cycles 
are absent in this case. The corona reaches a uniform energy state, which acts as an attractor 
in the diagram. 
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Fig. 7. — Photospheric azimuthal (rms) velocity amplitude with respect to time. The loop is 
subject to both nanoflare reconnection heating and Alfven wave heating. The amplitude of 
the Alfven waves increases with time becoming a non-negligible energy source in the second 
half of the simulation. Symbols indicate the same times as the corresponding same symbols 
in Fig. El 
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Fig. 8. — Phase diagram of the mean temperature and density in the corona in the case of a 
loop with both nanoflare reconnection heating and Alfven wave heating. Arrows indicate the 
time direction, and solid, dashed and dot-dashed curves (blue, green and red curves in the 
online version) denote limit cycles (the circle, triangle and square denote the start of these 
cycles, respectively). The dotted curve corresponds to the initial stage of the simulation. In 
this case the amplitude of the Alfven waves increases in time as indicated in Fig. [7J When 
the energy from the latter becomes non-negligible the corona reaches thermal equilibrium 
and correspondingly the cycles converge to a uniform state. 
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Table I. Falling velocities for coronal rain observed by Hinode/SOT 



Loop 


Apex 


Footpoint 


Free fall (i> =0) 


left 


50 - 70 


115 


115 


middle 


30 - 50 


80 - 140 


90 - 150 
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right 


60 - 70 


15 ->■ 105 


45 



Note. — Estimated velocities of coronal rain for 
the 3 loops of Fig. ??. "Left", "Middle" and "right" 
denote, respectively, the loops in the left, middle and 
right panels of that figure. The velocities are calcu- 
lated by estimating the slopes in the length - time 
diagrams. Speeds at the top part and bottom part of 
the loops are written in the "apex" and "footpoint" 
columns respectively. The "free fall" column is the 
velocity calculated assuming free fall and no initial 
speed. Positive values denote downward motions. In 
the loop of the middle panel an upward motion at a ve- 
locity of ~ 25 km s _1 is observed on the top part of the 
loop. In the loop of the right panel a deceleration to 
a velocity of ~ 15 km s _1 is clearly observed halfway 
down the loop. It then accelerates again to ~ 105 km 
s -1 . The free fall velocity for this case is calculated 
from the time of the start of this acceleration. 



